EPJ manuscript No. 

(will be inserted by the editor) 



Resonant charging of Xe clusters in Helium nanodroplets under 
intense laser fields 



o 



a 

O 

^ , 

(N 
> 

On 
O 
O 



Christian Peltz and Thomas Fennel 

Institute of Physics, University of Rostock, D- 18051 Rostock, Germany 
Received: date / Revised version: date 

Abstract. We theoretically investigate the impact of multiple plasmon resonances on the charging of Xe 
clusters embedded in He nanodroplets under intense pump-probe laser excitation (r = 25 fs, Io = 2.5 x 
10 14 W/cm 2 , A = 800 nm). Our molecular dynamics simulations on Xe3ogHeioooo and comparison to results 
for free Xe3og give clear evidence for selective resonance heating in the He shell and the Xe cluster, but 
no corresponding double hump feature in the final Xe charge spectra is found. Though the presence of the 
He shell substantially increases the maximum charge states, the pump-probe dynamics of the Xe spectra 
from the embedded system is similar to that of the free species. In strong contrast to that, the predicted 
electron spectra do show well- separated and pronounced features from highly efficient plasmon assisted 
electron acceleration for both resonances in the embedded clusters. A detailed analysis of the underlying 
ionization and recombination dynamics is presented and explains the apparent disaccord between the 
resonance features in the ion and electron spectra. 

PACS. 36.40.GK Plasma and collective effects in clusters - 36.40.WA Charged clusters - 52.25.Jm Ion- 
ization of plasmas - 52.50.Jm Plasma production and heating by laser beams (laser-foil, laser-cluster, etc.) 
- 52.65.Yy Molecular dynamics methods 
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1 Introduction 

Atomic clusters in intense laser fields provide a rich test- 
ing ground for exploring strong-field dynamics and col- 
lective phenomena in finite many particle systems [U[2]. 
Under intense near-infrared laser pulses (A ~ 800 nm), 
the transition of gas-phase clusters into well-isolated finite 
nanoplasmas results in high energy absorption, rapid clus- 
ter explosion, and the emission of highly charged ions [3 , 
[HH], fast electrons [HJEJEllH] •> an d energetic photons from 
the vacuum ultraviolet up to the x-ray range [l0l[TTJ[T2j 
fT3}fT4] . Because of the opportunity to study the underly- 
ing ionization, heating, and decay mechanisms in a nearly 
background- free environment, intense laser-cluster inter- 
actions are of high interest also for various other fields 
ranging from plasma physics to applied laser-matter re- 
search. 

Strong enhancements in the energy, yield, or charge 
states of emitted species [6,7,8,15 can be achieved with 
laser fields that induce resonant collective heating of nano- 
plasma electrons (Mie plasmon resonance). When present, 
the plasmon resonance is typically stronger than addi- 
tional higher-order nonlinear resonance effects [16 . How- 
ever, in rare-gas clusters where a dense nanoplasma is 
formed rapidly via tunnel ionization and subsequent elec- 
tron impact avalanche ionization, the high charge den- 
sity impedes resonant excitation in early stages of the in- 
teraction. The situation is similar for most simple metal 



clusters. Therefore, a certain cluster expansion is required 
for efficient resonant heating, making the temporal pulse 
structure to an efficient control knob for the dynamics. 
Control of the absorption [17] and the emission character- 
istics by the pulse structure has been demonstrated with 
various schemes, e.g. stretched pulses [4,6,8,15 , pump- 
probe techniques pT7 t fT8 ^ fT9 | [2Q] , and multi-parameter pulse 
shaping [5,21j. Typically, resonant plasmon excitations 
lead to the emission of laser aligned energetic electrons [6 , 
0[9] and highly charged atomic ions with high kinetic en- 
ergy [8,20j. Whereas the effect of resonance heating in ho- 
mogenous clusters is relatively well understood [22] , there 
is an ongoing debate on how the possibility of multiple 
plasmon resonances in core-shell systems affects the inter- 
action dynamics and the final emission spectra. Prominent 
specimens are clusters in He nanodroplets, as the pickup 
method with helium nanodroplets is well established in 
the strong- field cluster community [23p 4]. 

Recent theoretical work on Xe clusters in He droplets 
suggests that (i) nanoplasma generation is triggered by 
ionization of the embedded cluster and (ii) subsequent 
avalanche ionization rapidly ionizes the He nanomatrix [25j 
[26] . As the explosion of the He shell is fast, resonant con- 
ditions can be achieved at earlier times in the shell than 
in the embedded cluster. The above calculations show two 
separate resonance features in the energy absorption. One 
may expect that signatures of these multiple resonances 
should be imprinted also in the final emission spectra, e.g. 
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in pump probe experiments. So far, however, no indica- 
tions have been found in the published data to support the 
hypothesis that two different optimal delays - one tuned 
to a resonance of the rapidly exploding He shell and the 
other to the later resonance of the embedded cluster - 
can be used to enhance the emission of highly charged 
ions. The available data for embedded clusters appear to 
be qualitatively similar to experiments on free clusters and 
have thus often been interpreted by neglecting the detailed 
dynamics of the helium shell [19,23,27 . Nonetheless, it is 
by now well accepted, that a He shell has decisive impacts 
on the interaction process. A solid understanding of sig- 
natures from multiple resonances in core-shell systems is 
still lacking. 

In the current contribution we investigate this prob- 
lem theoretically within a pump-probe analysis using a 
classical molecular dynamics (MD) approach. The MD 
model accounts for tunnel and electron impact ionization 
of cluster and shell constituents under inclusion of local 
plasma field effects and considers electron-ion recombina- 
tion for the calculation of charge state distributions. As 
free and embedded model systems we consider Xe3og and 
XesogHeioooo under two linearly polarized 25 fs NIR pulses 
(800 nm) with intensity 2.5 x 10 14 W/cm 2 . Our systematic 
analysis of the spatially resolved energy absorption, the fi- 
nal ion charge state distributions, and the electron energy 
spectra as function of pulse delay reveals the following 
key results: (i) irrespective of the presence of the He shell, 
the final Xe charge states are substantially enhanced for 
a similar pulse delay and only one optimal delay is found 
in both cases. The amplitude of the enhancement is larger 
in the presence of the He shell, which can be attributed to 
more efficient inner ionization by the first pulse, i.e., prior 
to any resonance, and more efficient recombination due 
to stronger screening for nonresonant conditions, (ii) The 
spatially resolved energy absorption as well as the elec- 
tron emission spectra give clear evidence for the presence 
of two resonances with a pronounced double hump in the 
electron spectra as function of pulse delay. Our results 
thus indicate that the parallel measurement of electron 
and ion spectra offers a promising avenue to experimen- 
tally resolve multiple resonances in embedded clusters and 
their consequences for the ion spectra. 

The remainder of the text is structured as follows. In 
Sec. [2] we describe the theoretical approach and the nu- 
merical implementation of the MD model. Section [3] dis- 
cusses the key stages of the interaction dynamics based 
on a time-dependent analysis of representative simulation 
runs (Sec. l3.l() and the pump-probe spectra of the ions and 
electrons in relation to the spatially resolved absorption 
analysis (Sec. 13.2)) . Conclusions are drawn in Sec. HI 



2 Methods 

The explicit quantum mechanical treatment of clusters 
in intense laser fields is not feasible with current com- 
puter technology. We therefore employ classical molec- 
ular dynamics, which has been used with great success 




Fig. 1. Initial geometry of the investigated free (a) and embed- 
ded (b) clusters. In the embedded case, the relaxed icosahedral 
Xe309 is surrounded by a He shell in fee crystal structure. 



for describing laser-cluster interactions on a microscopic 
level, see e.g. Refs. [TJ[2] and references therein. The basic 
idea behind is that atomic ionization processes, which are 
strongly governed by quantum mechanics, are described 
statistically via appropriate rates, whereas the dynamics 
of the resulting ions and electrons can be treated clas- 
sically. A classical picture for plasma electrons and ions 
is in most cases justified for strong excitations resulting 
from intense laser irradiation. The laser-matter coupling 
is described in dipole approximation, which requires the 
clusters to be sufficiently small such that (i) the spatial 
variation of the laser field and (ii) laser propagation ef- 
fects in the target can be neglected. Both requirements 
are met in our scenarios. 



2.1 Molecular Dynamics Approach 

The following discussion summarizes methodical key as- 
pects of the dynamical model. Note that further details 
can be found elsewhere [26,33 . 

For the simulations we consider Xe3og in relaxed icosa- 
hedral geometry, see Fig. [TJi. A binary Lennard- Jones po- 
tential is employed for structural annealing [28] and switched 
off in the simulation of the laser-induced dynamics. The 
embedded XesogHeioooo cluster is constructed by cutting 
a sufficiently large spherical piece from an fee crystal with 
the density of liquid helium (p = 0.022 A~ 3 , [29J). Af- 
ter inserting the Xe3og core, He atoms with Xe neighbors 
closer than the equilibrium He-Xe distance of 4.15 A are 
removed [30], resulting in a compact core-shell compound, 
cf. Fig. [l]b. 

Inner ionization of atomic constituents is treated sta- 
tistically via rates for tunnel ionization (TI) and electron 
impact ionization (EII). Resulting ions and electrons are 
then propagated classically in the laser field and under 
the influence of mutual Coulomb interactions with other 
plasma particles. The particle trajectories are integrated 
by a velocity- Verlet propagator with 2.5 as timestep. The 
coupled equations of motion read 

rriiTi = ^eEi as - V Fl ^ (1) 



Christian Peltz, Thomas Fennel: Resonant charging of Xe clusters in Helium nanodroplets under intense laser fields 3 



where m^, ^e, and are the mass, charge, and position of 
the z-th particle and Ei as = e z Eo [f(t) + f(t — At)] cos(ut) 
describes the electric field of the two linearly polarized 
laser pulses with temporal delay At, respective peak am- 
plitude Eo = ^2Iq/cso, and 7kj=1.54eV photon energy 
(A =800 nm). Here Iq, c and eo denote the single pulse 
peak intensity, the vacuum speed of light, and the vac- 
uum permittivity. We employ a gaussian pulse envelope 
f(t) = exp(-21n2t 2 /r 2 ), where r is the pulse duration 
(FWHM). Note, that the single carrier wave results in 
constructive interference of the pulses for short delays. 
The pairwise Coulomb interaction Vij is described with a 
pseudopotential of the form 



Vi j (r ij ,q 1 ,q 2 ) 



QiQj 



erf 



(2) 



with the elementary charge e, the interparticle distance 
, their charge states qi and qj , and a numerical smooth- 
ing parameter s. The latter regularizes the Coulomb inter- 
action and offers a simple route to avoid classical recom- 
bination of electrons below the lowest possible quantum 
mechanical energy level. We use a constant value of s and 
choose it to be the minimal value, such that the binding 
energy of an electron to a g-fold charged ion is always 
lower than the q-th atomic ionization potential. In our 
case, the smoothing parameter is determined by Xe and 
has a value of s = 1.4 A. For the time consuming evalu- 
ation of the mutual particle-particle interactions we take 
advantage of massive parallel computation techniques. 

In each timestep, the probability for TI is evaluated 
from the Ammosov-Delone-Krainov rates [31] for each ion, 
employing the local effective electric field, i.e. the laser 
field plus the spatially smoothed local plasma field. The 
plasma field averaging avoids spurious TI events resulting 
from Rydberg electrons and is realized by evaluation of 
the pseudopotentials with a smoothing parameter similar 
to the ground-state Xe-Xe distance. A TI event results in 
a new plasma electron at the classical tunnel exit and an 
incremented charge state of the residual ion. 

Electron impact ionization is evaluated from effective 
Lotz cross sections [32], which have been modified to ac- 
count for local plasma field effects [33]. Considering an 
atomic ion in the charged cluster, neighboring ions and 
electronic screening by quasifree cluster electrons lead to 
an effective ionization threshold E* t = E n \ — Z\ e nv that is 
lowered with respect to the pure atomic value by the en- 
vironmental shift A env . Here E n \ denotes the bare atomic 
thresholds to remove an electron with principal and angu- 
lar quantum numbers n and I and E* t specifies the corre- 
sponding minimal energy required to lift the electron into 
the quasi-continuum within the cluster. The shift A env is 
evaluated directly from the plasma field in the MD simula- 
tion using the scheme developed in Ref. [33 . Therefore the 
local ionization barrier for a given ionic cell is assumed to 
be located at half distance to the next neighbor (rbar)- The 
shift Z\ env is then given by the electron potential produced 
from all charges (including the residual ion) at the point 
of the barrier (i*b a r) minus the electron potential in the en- 
vironmental field alone (residual ion removed) sampled at 




Fig. 2. Subdivision of the embedded cluster for the spatially 
resolved absorption analysis. Regions Q\ and i?2 contain the 
Xe core and the He shell, respectively (as indicated). The re- 
mainder of simulation space is denoted as the outside region 
J?3. Corresponding radii are updated during cluster expansion. 



the position of the ion. Whenever an electron penetrates 
a new ionic cell, the Lotz cross section cr E n is determined 
including all populated shells of the corresponding ion by 
employing the effective ionization thresholds E* t , the re- 
spective shell occupancies, and the kinetic energy of the 
impinging electron. Ionization is accepted, if the impact 
parameter b of the impinging electron fulfills the condi- 
tion b < ^fo^mfx. After a successful EII event the charge 
state of the corresponding atom or ion is increased and a 
new plasma electron is generated under the constraint of 
energy conservation. 



2.2 Spatially resolved energy absorption 

For the identification of resonance heating in different re- 
gions of the embedded cluster we perform a spatially re- 
solved analysis of the energy absorbed from the laser field. 
To this end, the embedded cluster is divided in three dif- 
ferent spherical regions, as indicated in Fig. [2j The core 
region Q\ contains the Xe cluster and the surrounding 
region i?2 contains the He shell. The rest is denoted as 
the outside region i?3. These regions are time-dependent 
and the corresponding radii are updated in each timestep. 
The instantaneous power absorption in region Qk can be 
expressed as function of the dipole velocity and the laser 
field by 

P Qk if) = e^rvE las (£), (3) 

rief2 k 



leading to the accumulated energy absorption 



P Qk (t')dt'. 



(4) 



The total absorption is then just given by the sum 

W ahs (t) = E W ^)- ( g ) 



k=1..3 



2.3 Determination of final charge spectra 

To determine the final ion charge spectra, recombination 
of quasifree electrons with cluster ions after the laser ex- 
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citation has to be taken into account. Note that the treat- 
ment of recombination has strong impacts on the charge 
distributions and is thus indispensable for predicting real- 
istic spectra [26) 133] . Two relevant mechanisms for electron- 
ion recombination during cluster expansion are radiative 
recombination and three-body recombination (TBR). As 
has been estimated previously [33 , radiative recombina- 
tion can be neglected due to low rates [53]. TBR, i.e. elec- 
tron capture after the collision of two quasifree electrons 
in the vicinity of an ion, proceeds mainly to high Ryd- 
berg states of the ion and can therefore be treated classi- 
cally. Hence, TBR is automatically included in reasonable 
approximation within the classical MD propagation. Be- 
sides the dominant contribution of TBR, even higher order 
collisional recombination processes (four-body, five-body, 
etc.) are accounted for because of the fully microscopic 
description of classical particle-particle correlations. The 
key advantage of the direct microscopic treatment of re- 
combination is that no approximations like quasi charge- 
neutrality or a thermal electron velocity distribution need 
to be used. In addition, also the local field effects from 
screening and potentials of neighboring ions are included. 

However, for a qualitative discussion of the recombina- 
tion dynamics, it is useful to recall basic approximations 
for the TBR rate from plasma physics. Considering a hy- 
drogen like ion in a high temperature plasma, the TBR 
rate for electron capture into a bound state with prin- 
cipal quantum number n is given by k n = an 4 TV 2 /T 2 , 
where N e and T are the number density and temperature 
of plasma electrons and a is a constant [35]. Summation 
over available n- levels up to the cutoff n c yields, in good 
approximation, a total TBR rate of 

k tot =anlN 2 JhT 2 . (6) 

With the thermal cutoff n c = n T = (i^i/Zc^T) 1 / 2 , where 
Ei is the binding energy for n = 1 and ks is the Boltz- 
mann constant, one finds the well-known fc t ot °c N% T -9 / 2 
TBR scaling law [36 . In contrast to that, when consid- 
ering a realistic cluster experiment, permanent recombi- 
nation is assumed only for electron capture in levels with 
n < n cr it, where n cr it indicates the critical level above 
which electrons can be re-ionized by the weak ion extrac- 
tion field of the detector (field ionization cutoff). 

When the highly excited cluster expands, the tempera- 
ture of quasifree electrons decreases rapidly via adiabatic 
expansion cooling (T ex N^ 3 ), leading to a steady de- 
crease of the TBR rate for fixed n c = n cr it : cf. Eq. (j6j). 
Hence, relevant n- levels with n < n cr n are efficiently pop- 
ulated only in early stages of expansion, while the total 
efficiency of this process is reduced for higher nanoplasma 
temperature, e.g., for resonant heating. Very high Ryd- 
berg levels (n > n cr i t ), on the other hand, are populated 
only at later stages, as high n- levels become available only 
after significant cooling and after a certain degree of ex- 
pansion (interatomic Coulomb barriers rise with cluster 
radius). These very high and long-living Rydberg levels 
are assumed to be re-ionized by the electric field of the 
ion detector or by other space charge fields. 



To measure the result of these complex recombination 
processes in the MD simulation, we apply the following 
simplified scheme. First, as the recombination dynamics 
is included in the MD, we simply propagate for a suffi- 
ciently long time after laser excitation. Electrons are then 
treated as recombined when localized to an ion (single par- 
ticle energy below the potential barrier of the respective 
ionic cell). The inner charge state of an ion is reduced by 
the number of localized electrons at the time t = t recom b. 
We use a sufficiently large value of t recom b, such that the 
resulting spectra are well-converged. The key assumption 
behind this simplified treatment is that n-levels leading to 
permanent recombination are populated only in a short 
time window after laser excitation. As has been checked 
carefully in a previous study [33] , recombination processes 
at longer time scales (for t ^> t reC omb) produce predomi- 
nantly very high lying Rydberg electrons that can be re- 
ionized easily by weak external fields. 



3 Results 

In the following we discuss the simulation results obtained 
with the above methods for pump-probe excitations of 
Xe309 and Xe3ogHeioooo by two 25 fs laser pulses at Iq = 
2.5 x 10 14 W/cm 2 . We begin with a time-dependent anal- 
ysis of typical runs to highlight the pivotal stages of the 
interaction processes in the free and embedded clusters 
and then proceed with the systematic inspection of ion 
and electron emission spectra. 



3.1 Time-dependent cluster response 

The evolution of selected key observables is displayed in 
Fig. [3] for three different pulse delays. Two of them are 
chosen such that resonant heating by the second pulse 
is induced in the He shell (left panels) or the Xe clus- 
ter (middle panels), respectively. The third case is a long 
delay, where the second pulse cannot excite either of the 
resonances (right panels). 

Focussing on the response to the pump pulse first, our 
analysis reveals the following picture: Irrespective of the 
presence of the He shell, ionization begins with tunnel ion- 
ization of Xe atoms, see the onsets of inner ionization in 
Fig. [3£l Laser heating of the first released electrons induces 
an impact ionization avalanche that quickly increases the 
inner Xe ionization levels. Up to the peak of the pump 
pulse, however, Xe inner ionization evolves very similar in 
the free and embedded systems, reflecting that the helium 
does not play a significant role in these very first stages. 
Subsequently, triggered by the plasma produced in the Xe 
core, also the He shell becomes rapidly ionized in a second 
EII avalanche to an average inner charge state of He ions 
of almost two. In turn, as a feedback effect of the shell 
activation in the second half of the first pulse, Xe inner 
ionization is enhanced via EII by almost a factor of two 
over the free cluster result (compare solid and dashed red 
curves in Fig. [3^). This effect can be attributed to the 
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Fig. 3. Calculated evolution of inner ionization, absorption, cluster and matrix radii (circles indicate estimated critical radii for 
resonant excitation), average Xe-core electron energy, and relative nanoplasma population in the Xe region (top to bottom) for 
free and embedded Xe309 (as indicated) under two 25 fs laser pulses (delays from left to right: At± — 100 fs, At<2 — 350 fs, Atz = 
800 fs) at intensity In = 2.5 x 10 14 W/cm 2 . Note the different scalings of the time axes. 



additional electrons produced in the He shell and the re- 
sulting stronger heating. The average kinetic energy in the 
Xe core region during the first pulse, which is crucial for 
EII of Xe, reaches an about three times higher peak value 
in the embedded case, see Fig.[3jl. After the first pulse, the 
total energy absorption of the embedded system is about 
20 times higher than for the free cluster, cf. dashed and 
solid black curves in Fig. [3]o. 

As a result of the ionization and heating induced by 
this first excitation step, the Xe cluster and the He shell (if 
present) begin to expand, see Fig. [3^. The shell explosion 
proceeds much more rapidly than that of the cluster due 
to the low atomic mass of He. Interestingly, the expansion 
speed of the Xe cluster changes only weakly in the presence 
of the He shell - see also Figs. [3h and[3hi, where a longer 
period can be inspected prior to the arrival of the second 
pulse. Only after several hundred fs a slightly slower Xe 
expansion is found in the embedded system. 

The He shell reaches its critical radius for resonant ex- 
citation after about t = 100 fs (see Fig.[3fc), which roughly 
coincidences with the arrival of the second pulse in the 
At = 100 fs case (left panels). The resonance heating of 
the He shell results in a total energy absorption that is 



by more than one magnitude higher than that from the 
first pulse (Fig. [3]o) . Note that the absorption in the core 
region (red curve in Fig. [3h) is still fairly low in this case. 
Nevertheless, though the Xe cluster has not reached its 
critical radius, absorption is also increased with respect 
to the first pulse (cf. the free cluster result in Fig. (3}d), 
but the enhancement is smaller than for the resonantly 
driven He in the embedded case. For the intermediate de- 
lay of At = 350 fs (middle panels), the Xe cluster (free 
and embedded) is near the critical radius at the arrival of 
the second pulse, see Fig. E3i. This is reflected in a much 
stronger absorption enhancement in the core region for the 
embedded system and in the total absorption of the free 
cluster (cf. Fig. [3^). Core absorption now exceeds 50% of 
the total energy absorption in the embedded case, reflect- 
ing efficient and selective resonance excitation of the Xe 
cluster. In contrast to that, only very weak heating by the 
second pulse is found for both cases (free and embedded) 
with the long delay of At = 800 fs (Fig. Eh), reflecting 
that no resonance can be excited. 

For all scenarios, Xe inner ionization stages are further 
increased by the second pulse, see Figs.|3h,|3f, and[3h, but 
the additional charging is always smaller than the gain 
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during the pump pulse. This behavior reflects that the 
EII cross sections rapidly decline with increasing charge 
states. Further, the additional inner ionization decreases 
with increasing pulse delay, which reflects that EII is most 
efficient if the Xe cluster is still dense. Inner ionization 
becomes less efficient as the Xe cluster expands, as both 
the flux of impinging electrons as well as the ionization 
threshold lowering decline. Final ion distributions are cal- 
culated by applying the recombination scheme described 
in Sec. 12.31 after a typical propagation time of 1 ps after 
the second pulse. Resulting average final Xe charge states 
are indicated in Figs.[3k,[3F, and [3k by the vertical arrows. 

The weakest charge state reduction by recombination 
is found with the intermediate pulse delay, where the Xe 
cluster is resonantly excited. Two effects are responsible 
for this behavior: First, as electrons in the Xe cluster are 
driven directly, the electron population in the Xe region is 
most strongly reduced by the probe pulse. This can be in- 
ferred from the evolution of the relative nanoplasma popu- 
lation in the core (see Figs. [3^, [3] and[3j)), which gives the 
ratio of electrons to total ion background charges in the Xe 
region. A lower population in the expansion phase reduces 
the recombination because of the density dependence of 
TBR. Second, the resonant heating strongly increases the 
electron temperature. The latter is measured by the mean 
kinetic electron energy in the Xe region, cf. Figs [3]i, [3j 
and [3fr. Comparing these average energies in the expand- 
ing cluster for a given Xe radius shows that the highest 
electron temperatures are achieved with At = 350 fs for 
the free as well as the embedded systems. This trend also 
suppresses recombination, because of the temperature de- 
pendence of the TBR rate, cf. Sec. 12.31 

For both the free and the embedded systems, a slightly 
larger effect of Xe recombination is found for the short de- 
lay of At — 100 fs, where the nanoplasma is less strongly 
depleted by the probe pulse and also the plasma temper- 
ature in the Xe region (for similar Xe radius) is slightly 
lower than for the At = 350 fs case. It should be em- 
phasized that even the strong resonant heating of the He 
shell with At = 100 fs is not sufficient to reverse the lat- 
ter trend. The strongest Xe charge state reduction by re- 
combination occurs with the long delay, where the high 
nanoplasma population and the low electron temperature 
in the core lead to efficient TBR. The particularly strong 
impact of recombination in the embedded cluster is at- 
tributed to the weak nanoplasma depletion. The latter 
effect reflects that the plasma in the surrounding He shell 
provides cold electrons to efficiently screen the inner Xe 
core. 



3.2 Pump-probe analysis 

In the next step we provide a more systematic analysis 
of the pump-probe effects by focussing on a smaller set of 
observables. Therefore the pulse delay has been scanned in 
a series of simulations and the resulting ion (upper panels) 
and electron spectra (lower panels) have been determined, 
see Fig. 0J The final emission spectra for the free (left 



panels) and embedded clusters (right panels) are shown 
in comparison to the energy absorption (middle panels). 

For both the free and embedded clusters, the final Xe 
ion distributions exhibit a pronounced pump-probe dy- 
namics (colored data in Figs. 0k and0Ji). The Xe charge 
states are maximized for similar optimal pulse delays of 
about At « 200 — 300 fs, producing substantially higher 
ionization levels than for much smaller or much larger 
pulse separations. The ion distributions are rather narrow 
with a typical spread between 3-5 charge states. 

Inspection of the results for free Xe3og reveals a pump- 
probe dynamics of the charge states that resembles the 
typical signatures of resonance enhancement. More pre- 
cisely, the optimal pulse delay for highest charging is in 
good agreement with the peak of the total energy absorp- 
tion (cf. Fig. 2k, andHJ}), reflecting plasmon enhanced ion- 
ization. The main reason for the pronounced pump-probe 
dynamics in the ion spectra is, however, not the varia- 
tion of inner ionization as function of pulse separation, 
but the strong effect of the delay on outer ionization and 
recombination. Comparison of the mean values for inner 
and final Xe ionization in Fig. 0k shows, that inner ioniza- 
tion is almost constant ((q^ e ) ~ 6.5) for all delays, with 
only a slight increase towards small delays (dashed line). 
The major effect for the enhancement at the optimal delay 
is the weak recombination at resonance, where the aver- 
age final Xe stages are close to inner ionization (compare 
dashed and solid curves). A clear marker for optimal res- 
onance excitation of the free cluster with At ~ 300 fs is 
the strongly enhanced emission of energetic electrons, see 
the electron spectra in Fig. 0k. Note that the feature at 
time zero can be attributed to constructive interference 
of the pulses whereas the asymptotic yield for long delays 
reflects the single pulse result. 

Focussing on the ion spectra from embedded clusters, 
a much stronger pump probe dynamics and still only one 
optimal delay are observed. At first glance, no clear effect 
from two resonances can be found in the ion spectrum, 
though the spatially resolved energy absorption in the He 
shell and in the Xe core indicate their selective excitations 
(cf. Fig. 0k). When compared to the free cluster, the max- 
imum final charge states at optimal delay are about 50% 
higher and ionization stages with long delays are lower. 
Inner ionization, though substantially increased over the 
free cluster case, is again only weakly dependent on the 
delay (when compared to the amplitude of the final charge 
states), but with a stronger enhancement at short delays. 
Note that the two small humps in the average Xe inner 
ionization at At = 100 and 300 fs may be interpreted as di- 
rect signatures from the shell and cluster resonances (cf. 
Fig. 0]i). Nevertheless, the stronger inner ionization in- 
crease at small delays is clearly ascribed to more effective 
heating of the helium that surrounds the still compact Xe 
core (density effect). 

Like in the free cluster case, most of the strong dy- 
namics in the final charge distribution is an effect of delay 
dependent recombination. The weakest charge state re- 
duction (cf. dashed and solid curve in Fig. 0]i) is found 
for a delay of At = 300 fs where the Xe core can be reso- 
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delay [fs] delay [fs] 

Fig. 4. Calculated final Xe charge spectra, energy absorption and electron energy spectra (top to bottom) for pure Xe309 (left) 
and Xe309Heioooo (right) under two 25 fs pulses (intensity: Io = 2.5 x 10 14 W/cm 2 ) as function of pulse delay. Average values of 
inner ionization (qa) (dashed) and final charge states (q^n) (solid) are given in panel (a) and (d). Solid curves in panels (e) and 
(f) represent the average kinetic energy (E^i n ) of emitted electrons. 



nantly excited. For long delays, recombination is particu- 
larly efficient in the embedded cluster, as the ionized He 
shell provides additional electrons for efficient screening 
of the Xe core. The maximum final Xe charge states are 
found for a delay between the optimal delay values for 
resonantly exciting the He shell or the Xe core, i.e. near 
At = 200 fs. This intermediate delay reflects a compro- 
mise between inner ionization, which is enhanced towards 
the He resonance, and weak recombination, which occurs 
close to the cluster resonance. 

Hence, besides a stronger pump probe dynamics and 
higher peak ionization levels, our model predicts no clear 
signature from multiple resonances in the final Xe spec- 
tra. In contrast to that, there is a pronounced double 
hump structure in the electron spectra, see Fig Sf. The 
two strong peaks in the maximum energy of emitted elec- 
trons near At = 150 and 370 fs can be directly assigned to 
plasmon enhanced electron acceleration [9 in the He shell 
and in the Xe core, respectively. Both coincide well with 
the absorption peaks in the respective regions, cf. Fig. 



4 Conclusions 

We have studied the ionization dynamics of Xe clusters 
embedded in He nanodroplets under strong-field pump- 
probe excitation to analyze the impact of multiple plas- 
mon resonances in the core-shell compound. 

The key findings are: (i) For the free and the embed- 
ded clusters we observe plasmon enhanced charging with 
pronounced enhancements of the final Xe charge states for 
similar optimal pump-probe delays, (ii) In the presence of 
the He shell, Xe ionization stages at optimal pulse delay 
are higher than for the free cluster, mostly due to more 
efficient inner ionization by the first pulse, (iii) The spa- 
tially resolved energy absorption analysis gives evidence 
that selective resonant excitation of the shell or the clus- 
ter is possible by a shorter or a longer delay, but without 
a corresponding double-hump feature in the ion spectra. 
Our analysis supports that the intermediate optimal de- 
lay for highest Xe charging reflects that both strong inner 
ionization and weak recombination are necessary to gen- 
erate high ionization stages, (iv) In contrast to that, the 
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two resonances do leave separate fingerprints in the elec- 
tron spectra, where highly efficient electron acceleration 
is found for the respective delays of the shell or core res- 
onance. 

These results thus indicate that a He resonance, even 
if present, may not be observed as a separate feature in 
the ion spectra. As direct cluster heating still plays a piv- 
otal role, a shell induced inner ionization enhancement 
only shifts the optimal delay for highest charging of clus- 
ter constituents. Our theory thus offers an explanation for 
the absence of corresponding features in the pump-probe 
experiments published so far. The strong features in our 
calculated electron spectra, however, suggest that exper- 
iments with parallel analysis on ion and electron spectra 
would offer a promising route to identify the presence and 
time-structure of multiple resonances for clarifying their 
link to the ion spectra. For such measurements very short 
pulses would be advantageous. 
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